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Generating multi-photon entangled states is a primary task for applications of quantum informa- 
tion processing. We investigate production of photon-triplet in a regime of light amplification in 
second-order nonlinear media under action of a pulsed laser beam. For this goal the process of cas- 
caded three-photon splitting in an optical cavity driven by a sequence of laser pulses with Gaussian 
time-dependent envelopes is investigated. Considering production of photon-triplet for short-time 
regime and in the cascaded three-wave collinear configuration we shortly analyze preparation of 
polarization-non-product states looking further applications of these results in the cascaded optical 
parametric oscillator. It is also demonststed the nonclassical characteritics of the photon-triplet in 
phase-space on the base of the Wigner function. Calculating the normalized third-order correla- 
tion functions below-and at the generation threshold of cascaded optical parametric oscillator, we 
demonstrate that in the pulsed regime, depending on the duration of pulses and the time-interval 
separations between them, the degree of three-photon-number correlation essentially exceed the 
analogous one for the case of continuous pumping. 



> 



PACS numbers: 03.65.Ud, 03.65.Yz 



I. INTRODUCTION 

Multiphoton entangled states have attracted a great interest in probing the foundations of quantum theory and 
constitute a powerful quantum resource with promising potential for various applications in quantum information 
technologies. Recently, experimental efforts in the direct production of multiphoton joint states, particularly, three- 
or four-photon states have paved a new stage for the study of multipartite entanglement [l[ . Indeed, the simultaneous 
generation of three photons is at the origin of intrinsic three-particle quantum properties such as Greenberger-Horne- 
Zeilinger (GHZ) -class and W-class quantum entanglement 0| - Up to now, several physical systems have been 
£ — , ■ proposed for the generation of photon triplet including third-order nonlinear medium [6( and cascaded spontaneous 
parametric down-conversion (PDC) 0, H| ■ Experimentally, three- photon down-conversion was studied in third-order 
nonlinear media § - 11 1 and also by using cascaded second-order nonlinear parametric processes [H)]. Direct gen- 



| eration of photon triplets using cascaded photon-pairs has been demonstrated in periodically poled lithium niobate 
crystals [l|. The distinction of three-photon GHZ and W states entangled in time and space has been also reported 
It was also shown that intracavity three-photon down-conversion can be effectively realized in cascaded 
optical parametric oscillator (OPO) [15j . This scheme that involves cascading second-order nonlinearities is based on 
the parametric processes of splitting and summing in which the frequencies between the pump and two subharmonics 
frequencies are in the ratio of 3 : 2 : 1. Experimental realization of cascaded OPO by using the dual-grid method 
of quasi-phase matching (QPM) has been done in Ref. [l6j]. Most recently, joint quantum states of three-photons 
with arbitrary spectral characteristics have been studied on the base of optical superlattises [ItJ for the cascaded 
configuration proposed in [To] ]. Two cascaded configurations have been considered in (l7j that lead to production of 
spontaneous photon triplet in cascaded PDC and generation of high intensity mode due to cascaded three-photon 
splitting in optical cavity. 

In this paper we continue the investigation of cascaded three-photon splitting in an optical cavity following the 
paper [17] . Our goal is twofold. In one part of the present paper we extend our previous results regarding three-photon 
splitting in optical cavity for an experimentally available scheme that is a cascaded parametric oscillator pumped by 
a sequence of Gaussian laser pulses (see, Sects. II and III). The other part of the paper is devoted to studies of 
quantum properties of "three-photon" mode. We discuss preparation of non-product states that are superposition of 
three-photon polarization states, however, without any consideration of cavity effects (see, Sec. IV). We also calculate 
the Wigner function of the subharmonic, i.e. "three-photon mode" showing the negativity in phase-space (see, Sec. 
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V). Our analysis includes calculation of third-order correlation function of photon numbers for various operational 
regimes of pulsed OPO(see, Sec. VI). It is known, that it is possible to control the behavior of quantum dissipative 
system by a train of pulses. In this paper, we use this approach for suppression of dissipation and cavity induced 
feedback in cascaded OPO that leads to increasing the level of three-photon-number correlation. 

II. PERIODICALLY PULSED CASCADED OPO: GENERATION THRESHOLD 

In this section we briefly describe the cascaded optical parametric oscillator (OPO) with a triply resonant optical 
ring cavity driven by a sequence of laser pulses. The semiclassical and quantum theories of this device for the 
monochromatic pump field were developed in Refs. [HI and here we only add some important details regarding 
cascaded OPO under laser pulses with Gaussian time-dependent envelopes. This cascaded configuration involves the 
fundamental mode driven by an external pump field at the central frequency ujq with an amplitude El and two 
subharmonic modes at the frequencies u)\ = ^ and u>2 = ^y 1 - Due to intracavity parametric type-I three- wave 
interactions pump field is converted to the subharmonics throughout two cascaded processes: loq — > uj\ + ui2 and 
L02 — > lo\ + w\. Subharmonic modes have the same plane polarizations and are all propagating in the same direction. 
The pump field consists from the sequence of Gaussian laser pulses with the amplitude 

E{t,z) = E L f{t)e- l ^> t - kLZ \ (1) 



/(t) = £e-(*-*<>— ) 2 / T2 , (2) 

n=0 

where T is the duration of pulses that are separated by time intervals r. 

The cascaded OPO is dissipative, because the modes suffer from losses due to partially transmission of light through 
the mirrors of the cavity. We consider below the case of high cavity losses for the pump mode (70 » 7,71 =72 =7), 
when the pump mode is eliminated adiabatically in non-depletion approximation. In this case, the effective interaction 
Hamiltonian in the rotating wave approximation reads as 

Hint = ihxiE f(t) (a£a% - aia 2 ) + ih\2 {al 2 a 2 - a\a^) , (3) 

where Eq = El/ Jo, a i (i=l,2) are the operators of the modes at the frequencies u>i = ^ and u>2 = ^f 1 and the 
coupling constants between the modes are expressed through the Fourier spectra of the second-order susceptibilities 
of nonlinear crystals of the length L 



Xi 



X2 



L 

dz X {2) (z)e lAkl{z)z , (4) 



n 



L 



dz X (2) (z)e lAk2 ^ z . (5) 







We assume collinear, one-dimensional on z quasi-phase-matching with the phase mismatch vectors 

Afci(z) = k L (uo,z) - ki(ujx,z) - k 2 (u>2,z), (6) 

Ak2(z) = k 2 (oj2,z) — 2ki(uJi,z) (7) 

analyzed in the details [TtJ • 

In this regime, the stochastic equations of motion for the complex c-number variables ot\ t 2 and flip corresponding 
to the operators a\, 2 and af 2 , have the following form 

^ = - 7iai + Eof{t)xifo + 2x 2 a 2 /?i + W ai (t), (8) 
= _ 7lA + E f(t) X ia 2 + 2x2/32^ + W Pl (t). (9) 
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The equations for a2,/?2 are obtained from (8), (9) by exchanging the subscripts (1) f± (2). Our derivation is based 
on the Ito stochastic calculus, and the nonzero stochastic correlators are: 

(W ai (t)W a2 (t))= Xl ^^S(t-t), (10) 
7o 

(W ai (t)W ai (t)) =2 X 2a 2 S(t -t), (11) 

(W Pl (t)W Pl (*')) = 2 X2 /3 2 <5(t - t), (12) 

(W aa (t)W aa {t))=2 X »a 1 8{t-t). (13) 

The Eqs. ([SI. (0) and the correlation functions modifies the analogous ones derived for OPO with monochromatic 
pumping [15j on the case of non-stationary pump field. 

In accordance with the cited paper, for the monochromatic driven cascaded OPO, zero-amplitude solutions a,\ — 
a 2 = of Eqs. ©, © with f(t) = 1 are stable, if El < E t h = ^v / 7i72, while steady-state photon numbers n\ and 

n 2 display histeresis-cycle behavior in a small domain < El/E^ < 1. Thus, remarkable feature of OPO under 
monochromatic pump is comparatively low generation threshold that depends from the second-order susceptibility in 
comparison with the scheme of direct intracavity three-photon down-conversion, where the pump power threshold is 
determined by the third-order susceptibility Q. 

Below, we derive the threshold value for OPO driven by trains of Gaussian pulses. The analysis of stochastic 
equations shows that similar to the standard OPO with monochromatic pump field amplitude, the periodically pulsed 
OPO also exhibits threshold behavior, which is easily described through the period averaged pump field amplitude 
f(t) = — L f(t)dt. We demonstrate this statement analyzing the stability of zero-amplitude solutions of ot\ — a 2 — 
of Eqs. (J8J), ([9]) for both modes below threshold and for the case when the decay rates of subharmonics are equal 
one to the order, 71 = 72 = 7. To check the stability we turn to the linearized on the small deviations 5cti,6f3i the 
equations in the semicalssical approach without noise terms. These equations can be rewriten in the following form: 

±5X ± = (- 1 ±^ Xl f(t))SX ± , (14) 
at 70 

±5Y ± = (- 7 ±^xif(t))6Y±, (15) 
at 70 

through the quadrature field variables defined as Sa± = 8ot\ ± 5a 2 and SX± = h(Sa± + 6a±), 5Y± = ^(Sa± — 5a±). 
In these variables the time evolution has the simple form 

6X±(t) = A±(t,t )5X ± (t ), (16) 
SY ± (t) = A T (t,t )5Y ± (t ), (17) 



A±(Mo) = ex P {±^± f f(t)dt - 7 (t-t )). (18) 
7o Jo 

Analyzing semiclassical equations and operational regimes we choose the switching time in infinity, i.e. to — > —00, 
and add in (J2j) terms with negative n. In this case, the function f(t) is periodic on time f(t + r) = f(t) and the 
analysis is simplified. Since the function f(t) is periodic on time, we can obtain the general formula J t f(t)dt = 
f(t)(t - t ) + F(t) - F(t ), where F(t) is a periodic function, F(t + r) = F(t). Therefore, we see from ([IB]). ([17). 
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(fT8| that the solution a; = below-threshold is stable if El < ~^-f{t). It is easy to check also that due to noted 
periodicity of the amplitude the following formula takes place (see, for example |18| ) 

/ dt J2 e- ( *-" r)2/T2 = / dte- t2 / T \ (19) 

"'O n=-oo 



This formula allows us to calculate the averaged value of the amplitude f(t). On the whole, we arrive to the result 
that for the case of Gaussian pulses above threshold regime is realized if 

E L >E^=-^=E th = -^=^l. (20) 

The important peculiarity of the system proposed is that the threshold value E t h depends on the coupling constant 
Xi which is related to the second-order susceptibility as well as depends on the characteristics of laser pulses. 



III. NUMERICAL SIMULATION OF DISSIPATION AND DECOHERENCE 



The cascaded OPO is dissipative, because the modes suffer from losses due to partial transmission of light through 
the mirrors of the cavity and due to quantum fluctuations. We analyse dissipative and decoherence effects on the 
base of master equation for the density operator of the cavity modes in the Lindblad form 

d 1 

-qT = -^[H in t,p] + ^2 H faipaf - ajaip - pa+a;) . (21) 

' ' ' ' i=l,2 

We calculate the quantities of interest (the photon number distributions, Wigner functions, etc.) mainly for the 
subharmonic mode (1) by using the reduced density operator p±(t) which is constructed from the density operator 
p(t) of both modes by tracing over the mode (2), pi(t) = Tr2{p). We analyze the master equation numerically using 
quantum state diffusion method (QSD) [l9j], [20]. According to this method, the reduced density operator is calculated 
as the ensemble mean 



1 N 

p(t) = M(\^(t))(^(t)\) = lim J2\Mt))mt)\ (22) 

N— >oo iv * — ' 



over the stochastic pure states \ip^(t)) describing evolution along a quantum trajectory. The stochastic equation for 
the state \ip^(t)) involves Hamiltonian described by Eq. ([3]) and the Lindblad operators described by noise terms in 
the master equation (|2"TT) . We calculate the density operator using an expansion of the state vector |^) in a truncated 
basis of Fock's number states of modes of the subharmonics (1) and (2) 

IM«))=E a ™(*)l^>il^>2. (23) 

n 

Details of analogous calculations for an anharmonic oscillator in time- modulated field can be found in [2l|). The 
numerical simulations are performed in the truncated Fock basis of the subharmonic modes that are limited by 
500 photons. This approximation is valid for the case of strong nonlinear couplings \i and %2 with respect to the 
dissipation parameters. 



IV. PRODUCTION OF POLARIZATION, NON-PRODUCT STATES OF PHOTON TRIPLET IN 

COLLINEAR CONFIGURATION 



Three-photon correlations allow the creation of tripartite entangled states such as the GHZ state. For cascaded 
SPDC in noncollinear configuration spatially-polarization GHZ states has been considered in fl7l ]. Below we apply 
the results obtained for preparation of three-photon polarization states in collinear configuration of interacting waves. 
Recently, a simple but highly efficient source of polarization-entangled photon pairs at nondegenerate wavelengths and 
in collinear configuration has been demonstrated [22[ • We consider the production of polarization-entangled photon 
triplet. It is possible for the case when cascaded processes involve polarized photons. 

Thus, we modify the above results considering three-wave interaction with the indexes of polarization states. 
Looking further applications of above results for intracavity three-photon down-conversion in collinear configuration 
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of cascading processes we concentrate on consideration of non-product states that are entangled only on polarization 
degree of freedom but not on spatial variables. Thus, including into consideration the polarization states of the 
photons we assume that the type-II process wo ->wi+W2 create the pair of photons with a vertical V, (a% (V), (a2(V)) 
and horizontal H, (ax(H), ((12(H)) polarizations in collinear configuration. If the pump field is oriented at 45° to the 
horizontal and vertical axes two processes ujq — > u)\ (V) +W2 (H) and luq — > lj\ (H) +LU2 (V) take place in the first crystal. 
The next process is considered as the type-I parametric process. In type-I conversion, photon pairs are created with 
the same polarization state, but ortogonal to the input mode. Therefore, the second, type-I crystal is arranged in the 
manner that the following process: 0J2(H) —> oj\(V) + uj\(V) and ui2(V) —> ui\(H) + uj\(H) should be realized. 

For simplicity, we restrict our attention considering frequency-uncorrelated three-photon states and assume that 
the process under photons with (V) and (H) polarizations are described by the equal coupling constants. Thus, we 
assume that photon pairs in three- wave processes : ujo — > ui\ + ui2 , ^2 — > u l + uj 1 have correlations on the polarization, 
but not on the spectral lines. In analogy with Eq.(3), we model the sum of the corresponding parametric interactions 
by the following effective Hamiltonian 

H = H l +H 2l (24) 
Hi = ih X E Q f(t)(a+(V)a+{H) + a+(H)a+(V)) + h.c, (25) 
H 2 = ihk[{a 2 (V){a+(H)) 2 + a 2 (H)(a+ (V)) 2 )] + h.c. (26) 

Here, di(V) and 02(F) are the annihilation operators of modes (1) and (2) at vertical polarizations, while the operators 
a\(H) , 0.2(H) corresponds to the horizontal polarized photons of the frequencies uj\ and respectively. 

We focus on analysing the generation of non-product states for one-passing configuration of cascaded parametric 
spontaneous processes without consideration of cavity dynamics and feedback effects. This approach is valid for short 
interaction time intervals much shorter than the characteristics relaxation time. In this case time-evolution of the 
vector state of the system is described by the second-order term of the perturbation theory. Choosing the initial 
state as a vacuum state \ipi n ) = |0) ai (m|0) a2 (m|0) ai (y)|0) a2 (y) for all modes we derive the final state during time 
evaluation in the following form 

\m) = (--) t 2 H 2 H 1 \^) in = 

X~kE (\V)\V)\V) + \H)\H)\H))\0) a2{H) \0) a2{v)l (27) 

where \ = XP an d k — kt are the coupling constants and the states \V) = af(V)\0) ai(v) , \H) = (H)\0) ai(H) present 
the vertical and horizontal polarization states of photons at the frequency u)\ = Thus, we demonstrate that in this 
collinear, one-dimensional cascaded scheme triple photons can constitutes the polarization entangled (non-product) 
states of light. It should be noted that under cavity feedback effects the non-product quantum state cannot be 
described by this simple expression and we should include the higher-order terms of the perturbation expansion into 
consideration. 



V. PHOTON TRIPLET IN PHASE-SPACE: WIGNER FUNCTIONS AND PHOTON NUMBER 

DISTRIBUTIONS IN THE PULSED REGIME 

Quantum interference signature of three-photon states inphase-space has been demonstrated for the direct three- 
photon down-conversion in third-order nonlinear medium [y, Q as well as in the cascaded scheme [ItJ for the case 
of monochromatic pumping. We demonstrate now this effect for the pulsed regime of cascaded OPO. We illustrate 
these effects numerically on the base of the master equation, however, in the regimes when the dissipation in the 
cavity is unessential and the dynamic of modes is almost unitary. For the cavity configuration presented, the validity 
of such approximation is guaranteed by consideration of short interaction time for which the duration of pulses are 
much shorter than the characteristics relaxation time, l/(xi-Eo), I/X2 << t << 1/71,2 , provided that the nonlinear 
coupling constants exceed the dumping rates for the modes. 

Below, we present the results on the photon-number distributions and the Wigner functions for three-photon mode. 
The photon number distribution for oj\ = ^ mode is calculated as the diagonal element P\(n) — (n\p\n) of photonic 
Fock states while calculations of the Wigner function for the this mode are performed by using its standard formula 
in a Fock space: 

W 1 ( P ,0) =Y,Pl,rnnW mn {p,0). (28) 

771,71 
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FIG. 1: The photon number distributions (a, d), the Wigner functions (b, e) and its contour plots (c, f) for ui-mode and for 
different time intervals: t = to + 6.2 * 10 _3 7 _1 (a, b, c); t — to + 6.2 * 10 _2 7 _1 (d, e, f). The other parameters are: X1/71 = 200, 
X2/71 = 100, 72 = 71 = 7,T 7 = 1* 10~ 3 ,T7 = 6* 10~ 3 . 



Here, p, 9 are the polar coordinates in the complex phase space which is determined by the position and the momentum 
of quadratures x = (a + a + )/v / 2, y = (a — a + )/^/2i, respectively, while the coefficients W mn (p,0) are the Fourier 
transforms of the matrix elements of the Wigner characteristic function. 

Three-photon structure of the mode is shown on the photon- number distributions in Figs. 1 (a, d). As we see, for 
short time- intervals, the most probable values of photon numbers are separated by three photons. In Figs.l (b, e) we 
present results on the Wigner functions that clearly exhibit phase-space quantum interfringes. These results describe 
the case of cascaded OPO under two consecutive pulses with the duration T7 = 1* 10~ 3 separated by the time-interval 
tj = 6 * 1Q~ 3 . Fig. 1(b) shows the Wigner function evolved for a time interval t = to + 6.2 * 10 _3 7 _1 that corresponds 
to maximal photon number n max of the first pulse; Fig.l (e) shows the Wigner function at t = to + 6.2 * 10 — 2 7 — 1 
corresponding to n max of the second pulse. The Wigner function show three phase components with an interference 
pattern in the regions between them. We show the regions of quantum interference in the contour plots (see, Figs.l 
(c, f)) depicting negative regions of the interference terms in black. Note that threefold symmetry of the Wigner 
function and interference pattern has been demonstrated for the direct three-photon down-conversion in x^ 3 ' media 
0, Q. However, we note that the results presented here for the pulsed cascaded configuration are also different in 
details from the analogous calculation of the Wigner function for the case of monochromatic pump field (l7| . 



VI. PHOTON-NUMBER CORRELATION IN THE PULSED REGIME 



The experimental verification of time-dependent correlation between photons in triplet has been demonstrated for 
one-passing configuration of cascaded SPDC [l|. Considering production of photon triplet in a cavity, it seems that 
the correlation between photons can be evidently displayed for short interaction time intervals much shorter than 
the relaxation time. Nevertheless, the three-photon number correlation exceeding the coherent level, (that means 
the normalized third-order correlation function </ 3 ' > 1), has been demonstrated for the cascaded OPO driven by 
monochromatic pumping, in over transient regime for modes generated below the threshold (l7| . This effect decreases 
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if the system moves to the range of the generation threshold. At the threshold, the typical value for the normalized 
third-order correlation function for zero delay-time, — 1.2 has been obtained. 

Note, that effects of two-photon correlations for ordinary OPO and NOPO have been a subject of intense research 
efforts for last years. These problems have been investigated in details in a linear approximation on quantum fluc- 
tuations as well as in the framework of exact quantum theory of intracavity parametric generation with allowance 
for quantum noise of arbitrary intensity. In this way, the critical behavior of the second-order correlation functions 
which describe photon correlation effects has been found analytically on the base of the Fokker-Planck equation for 
the density matrix in the threshold region of generation [23| . 

In this Secion, we demonstrate the new regimes of strong three-photon correlation for the pulsed cascaded OPO. We 
concentrate on the numerical simulation of both the mean photon number of subharmonics as well as the third-order 
correlation function. Let us now discuss photon-number correlation in the time domain, considering output twin light 
beams from the pulsed cascaded OPO on the base of normalized third-order photon-numbers correlation function g^ 
for the mode (1) 

g^(t)=Tr(at 3 al Pl (t))/n 3 (t). (29) 

Here, n(t) = Tr{a J [a\p\ (t)) is the mean photon number. Considering three-photon number correlation for the intensive 
cavity mode in the presence of dissipation and cavity induced feedback, we control quantum dynamics of dissipative 
systems by the train of pulses. Here, we use this approach for controlling quantum statistics of mode, particularly 
for increasing the level of three-photon- number correlation. We analyze the cases in which T7 < 1 and T7 > 1, for 
over transient time-intervals, t 3> 7 -1 , considering the operational regimes below-and at the generation threshold. 
Typical results for the mean photon numbers and the correlation function g^ for two different parameters of the 
Gaussian pulses : T = 7~ 1 ,r = 107 _1 and T = 7~ 1 ,r = 57 -1 are presented in Figs. ([2} and As we see from 




FIG. 2: The mean photon number (a) and the third-order correlation function (b) versus £7 at the threshold EL/Eth = 1, 
curves (1) and below the threshold Ei/E t h = 0.5, curves (2). The parameters are: X1/7 — 0.2, X2/7 = 0.1, 72 = 71 = 7, 

TJ = 10, T7 = 1. 

these figures, the time-dependence of these quantities repeat the periodicity of the pump laser over transient time 
intervals. We also conclude that the maxima of the three-photon correlation are realized for the definite time intervals 
for which the mean photon number of mode (1) is in the ranges between its maxima and minima. As shown our 
calculations, such strong photon correlations take place in nonstationary regime of cascaded OPO when duration of 
pulses is close to a characteristic dissipation time. In Figs. ([2]) we analyze the mean photon number (Fig. [2ja)) and 
the correlation function (Figj2][b)) in the operational regimes of OPO below- and at the generation threshold. In the 
regime below the generation threshold, at E/Eth = 0.5 the correlation function display two-peak structure (see, curve 
(2)). The lower peak corresponds to the time intervals for which the mean photon number is between its minimal and 
maximal values. More correctly, for the period trains of intervals t = <o + Jn.ll. 67 , (m — 1, 2, ...), the mean photon 
number n = 50 and the correlation function g( 3 ' = 8.9. The other peaks correspond to the minimal values of the 
mean photon number. At the threshold (see, curves (1)) the effect of photon correlation is decreased, although the 
level of correlation exceeds the coherent level, g^ > 1, particularly, we get g^ = 3.4 for the mean photon number 



8 



16-, 




FIG. 3: Mean photon number (a) and the correlation function (b) for continuously pumped 0P0 at the threshold El/E^ = 1- 
The parameters are: xi/7 = 0-2, X2/7 = 0.1, 72 = 71 = 7. 

n=251. Thus, we found a remarkable result that the degree of three-photon number correlation for the pulsed regime 
of OPO surpass the analogous result for OPO with continuous pumping for the same mean photon numbers. Indeed, 
this conclusion is illustrated in Figs. ((2} and ([3]), where the comparison of the results on pulsed regime with the 
calculations based on the Hamiltonian ([3]) with / = 1 is done. Note, that the ideal limit of continuous pumping is 
realized if T — > 0, t — > 00 for the case of infinity numbers of pulses. We present on the Figs. ([3]) the results for OPO 
with continuous pumping at the threshold El = E t h, in which the mean photon number n = 52 in the steady state 
regime (curve (2), Fig. EJa)) approximately equals to the maximal value of the photon number in the pulsed regime. 
However, as we see, the level of the maximal correlation, — 9 in this case exceeds the analogous one for the case 
of continuous pumping, g^ — 1.2. 

It is natural to explain such improvement of three-photon correlation by control the behavior of a quantum system 
by an external time-dependent force. The presence of these effects in the cascaded OPO, particularly, can been seen 
from the noise-correlation functions. Indeed, the equation (|10p describes a multiplicative noise-term, where the level 
of noise is determined by the amplitude of pulsed driving field Ei,f(t) leading to the control of dissipation. In this 
spirit, we emphasize that the idea of controlling the dynamics of a quantum system in the presence of dissipation and 
decoherence by an external periodic driving was exploited by many authors (see, for example, [24] and the references 
therein). In one of the standard techniques control of the optical quantum system is achieved through the application 
of suitable tailored, synchronized laser pulses 25]. In this way, it is interesting to analyze the three-photon correlation 
function in dependence of the time-separation between pulses, i.e. for the other parameters of driving pulses in 
additional to the parameters considered in Figs. [5J The results are presented in Figs. 0] in the regime below the 
threshold where strong three-photon correlations are realized. As we see, in this case the correlation function reach 
to f^ 3 ' = 75 for n — 1.1 (curves (1)) for time-intervals between pulses, r = 107 _1 . Decreasing of the time-separation 
between pulses leads to decreasing of the correlation function (see, curves (2), where r = 57 _1 ). 

VII. CONCLUSION 

In conclusion, we have studied quantum properties of photon triplet cardinally different from those of twin photons. 
Because photon-triplet originates from a single laser photon, the quantum correlations take place between all three 
photons allowing the creation of entangled, non-product states. The production of photon-triplets in the presence of 
stimulation radiative processes, cavity feed-back effects and dissipation have been investigated. We have demonstrated 
the possibility to create polarization, non-product states of photon-triplet for one-passing, collinear configuration 
of cascaded parametric spontaneous processes. We have also illustrated three-photon structure of sub-harmonic 
mode on the base of both the photon-number distribution and the Wigner function. We have demonstrated the 
operational regimes depending on the durations of pulses and the intervals between them that guarantees strong three- 
photon-number correlations. This effect of strong correlation takes place for the definite time intervals corresponding 
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FIG. 4: The correlation functions (a) and the mean photon numbers (b) at EL/Eth = 0.2 for time-separation between pulses: 
T7 = 10, curves (1) ; T7 = 5, curves (2). The parameters are: X1/7 = 0.2, X2/7 = 0.1, 72 = 71 = 7, T7 = 1. 

to generation of high intensity "three-photon mode" in over transient regime and for wide ranges of the system 
parameters. We hope that these results could be of interest in areas of quantum communications and photonic 
quantum computing. 
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